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^—1 I Abstract 

CN ' Following the recent recognition of a positive value for the vacuum energy density 

^*p. . and the realization that a simple Kantowski-Sachs model might fit the classical tests 

l/^ I of cosmology, we study the qualitative behavior of three anisotropic and homogeneous 

^^ ■ models: Kantowski-Sachs, Bianchi type-I and Bianchi type-Ill universes, with dust 

^— V . and a cosmological constant, in order to find out which are physically permitted. We 

^^ I find that these models undergo isotropization up to the point that the observations 

^^ ' will not be able to distinguish between them and the standard model, except for the 

O . Kantowski-Sachs model {^ko < 0) and for the Bianchi type-Ill {^ko > 0) with f^Ao 

I [ smaller than some critical value ^Aj^j ■ Even if one imposes that the Universe should be 

trv| nearly isotropic since the last scattering epoch {z ~ 1000), meaning that the Universe 

H" ■ should have approximately the same Hubble parameter in all directions (considering 

• 1— I , the COBE 4- Year data), there is still a large range for the matter density parameter 

K> I compatible with Kantowsky-Sachs and Bianchi type-Ill if |ilo + ^Ao ~ 1| ^ ^^ foi^ ^ very 

C^ ■ small 5 . The Bianchi type-I model becomes exactly isotropic owing to our restrictions 
and we have Oq + ^a^ = 1 in this case. Of course, all these models approach locally 
an exponential expanding state provided the cosmological constant J7a > ^Am • 

1 Introduction 

Over the last five years, the issue of whether or not there is a nonzero value for the vacuum 
energy density, or cosmological constant, has been raised quite often. Indeed, the possibility 
of a nonzero cosmological constant A has been entertained several times in the past for 
theoretical and observational reasons (the early history of A as a parameter in General 
Relativity has been reviewed by [|ll, @, and 0). Recent supernovae results [Q, [|1 strongly 
support a positive and possibly quite large cosmological constant. Even taking the Hubble 
constant to be in the range 60-75 km/s/Mpc it is possible to show 0] that the standard 
model of fiat space with vanishing cosmological constant is ruled out. In a very nice review 
it is argued that postulating an ^A-dominated model seems to solve a lot of problems at 
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once. And again, in a quite recent review on the physics and cosmology of the cosmological 
constant, it is added that "recent years have provided the best evidence yet that this elusive 
quantity does play an important dynamical role in the universe" [§]. 

On the other hand, if the classical tests of cosmology are applied to a simple Kantowski- 
Sachs metric and the results compared with those obtained for the standard model, the 
observations will not be able to distinguish between these models if the Hubble parameters 
along the orthogonal directions are assumed to be approximately equal 0. Indeed, as JlO 



points out, the number of cosmological solutions which demonstrate exact isotropy well after 
the big bang origin of the Universe is a small fraction of the set of allowable solutions of 
the cosmological equations. It is therefore prudent to take seriously the possibility that the 
Universe is expanding anisotropically. Note also that some shear free anisotropic models 



display a FLRW-like behavior, as it is shown in |11 



2 The global behavior of the A 7^ models 

Taking all this into consideration, we discuss the behavior of some homogeneous but anisotropic 
models with axial symmetry, filled with a pressureless perfect fluid (dust) and a non vanishing 
cosmological constant. For this, we will restrict our study to the the metric forms 

ds^ = -c^df + a\t)dr^ + h\t) [de^ + fl{e)d<p'') , (1) 

with the two scale factors a{t) and 6(t); k is the curvature index of the 2-dimensional surface 
dO"^ + fki^)d(t>^ and can take the values +1,0,-1 implying fk{d) equal sin(6'), 6, sinh(6'), 
respectively, giving the following three different metrics: Kantowski-Sachs, Bianchi type-I, 
and Bianchi type-Ill [|12| , |1^ . 



Einstein equations for the metric (|T]), for which the matter content is in the form of a 



perfect fluid and a cosmological term. A, are then as follows |T^, |]T3|: 



^ a h IP' kc^ ^ ^ , 2 /r,\ 

2 - I + l2 + i;r = S^G'p + Kc\ 2 

a 0^ 0^ 

^\4^'^--^^GU^c^ (3) 

- + 7 + - 7 = -SttG- + Ac^ 4 

a b a b c 

where p is the matter density and p is the (isotropic) pressure of the fluid. Here G is the 
Newton's gravitational constant and c is the speed of light. If we consider a vanishing 
pressure {p = 0), which is compatible with the present conditions for the Universe, the last 
two equations take the form 

2^ + (^ + T^ = ^^- (^' 

« + L«^ = Ac^ (6) 

a b a b 

and equation (^ can easily be integrated to give 

fe2 _ Ml ^ Ax^ A 2 



where Mi is a constant of integration. 

The Hubble parameters corresponding to the scale factors a{t) and b{t) are defined by 

Ha = a/a and Ht = b/b. 

Using them to introduce the following dimensionless parameters, in analogy with which 
it is usually done in the Friedmann-Lemaitre-Robertson- Walker (FLRW) universes, let us 

define 

Ml 



b'H, 



2 = ^M, 



^^n, (9) 



bm, 



and 

3i?j-"- (10) 

The conservation equation (^ can now be rewritten as 

VtM + fifc + l^A = 1. (11) 

Now defining the dimensionless variable y = b/bo where 6o = ^(^o) is the angular scale 
factor for the present age of the Universe, and using equation ([TT| ) (taken for t = to), one 
may rewrite equation (0) as 



y = ±H,JnMo{^ - 1) + fiAo(l/' - 1) + 1, (12) 

where the density parameters defined previously and Hf, with the zero subscript, denote 
as before these quantities at the present time to- Iii this way, the number of independent 
parameters have been reduced. Substituting equation (^ into equation (||) gives 

. ^ M, - Ml f + lAcV 
2^Mib-kc%^ + ^c%^' 

where Mp is a constant proportional to the matter in the Universe, 

Mp = SnGpab"^. (14) 

Using the procedure above, equation (|T3D can be rewritten in the form 

np-nM + 2nA = 2^, (15) 

where 



^^ - ^- ^''^ 



__M„ 

From equation (0) one may define a matter density parameter. For this, we introduce the 
notion of mean Hubble factor H such that 3H = Ha + 2Hf,. Also, for these models, the shear 
scalar a |T^ is given by 

a=^iHa-Hb) (17) 

3 



Thus, equation (^ may be rewritten [12 



as 






kc\ 



As in FLRW we call critical matter density pc when A; = and A = 



Pc 



3H^ 



a 



SnG 



The matter density is generally defined as i7 = p/pc, then 

SnGp _ 8nGp 



n 



3/72-^2 2HaHi, + Hf 



(19) 



(20) 



just like in FLRW models, and such that f2 = 1 when A; = and A = 0, and which is related 

to Qp by 

Q = -^. (21) 



1 + 2 



Ha- 
Hh 



Although Qm is not the matter density parameter, it performs the same important role. We 
emphasize the fact that if for one particular time Ha = H^ and Vt^^ = 1, then, by equations 
dnl), (|1|) and (HD, 3fi = fip = ^m = -^k] and if < ^a < 1 and ^m = 1, then, -fi^ = ^a 
and fi ~ 1. 

Introducing another dimensionless variable x = a/aQ, equation ([T3|) takes the form 



X 



H, 



^(l-f)+fiAo(-l + V) + i^ 



bo 



bo 



WfiMo(^-l)+fiAo(y'-l) + l 



f22) 



and its number of independent parameters was also reduced, now at the expense of equation 
([T5| ) taken for the present time t = to- 

Now, we want to find the time dependence of h{t) in a qualitative way, starting from 
equation (|1^). Since the model universe will be defined only where ij'^ > 0, as was previously 
done by |n| for FLRW models, the problem is reduced to find out the zeros of y, with y y^ 0. 

There are two Q\ values which characterize two zones of distinct behavior for the scale 
factor b. Starting with condition y = one may obtain 



^Ao 



{^mq -l)y-^ 



Mo 



y" -y 



(23) 



If we consider I^Aq = QAoiy)^ as a function of y, then this function presents a relative 
minimum and a maximum, that we will denote by ^a^ and ^Am, respectively. The relative 
minimum depends on Qmo ^^ the following way fl^: For Qmo < 1/2 we have 



Q: 



3fi 



A/o 



\ 



"Afo 



-1 



Q 



A/o 



n 



Mo 



1/3 



+ 



^(Omo - l)V^i/o - 1 + (1 - ^Mo)/^Mo 



1/3 



(n 



Mo 



(24) 



for Qmo > 1/2 the expression is 



^A. = -SQmo cos f ^^ ) - (f^A/„ - 1). (25) 



The relative niaximum is done by 



^Am = -S^A/o cos I — - — 1 - {Qmo - 1), (26) 

where 9 = cos~^ ( -^ — j. These expressions are hmiting zones of the (^Aq, fijv/o) plane, 

where y = has three or one solutions (for details see [|14|)- The ^Am expression is also 
defined for Qmo > 1/2? but it has the meaning of a maximum only for Qmq > 1- The i^A,, 
less or equal to ^Am imposes the recollapse of scale factor b, while greater values produces 
inflexional behaviors for b. The Qaq values greater or equal to ^a^ are physically "forbidden" 
because they don't reproduce the present Universe (see [|14|)- Obviously, ^a^ < ^Ac always. 
Although we are considering anisotropic models, the equation ([T2|) for y as a function of 
Qmo is mathematically the same as equation (2) obtained by |]1^ for the homogeneous and 
isotropic FLRW models. From equations (0) and (^) one obtains the differential equation 



n. 



(27) 



dx ^a-f)+^Ao(-i+^^^) + ^,, 

dy ^Mo ( 1 - y) + ^Ao (y^ -y)+y 

This equation automatically complies with the two conservation equations ( p^lD and (|I5|) 
evaluated at to- There are some particular values of the parameters {Qmo,^Ao) for which 
this equation has exact solutions. However, for the majority of the values of the parameters, 
the solution has only been obtained by numerical integration. 

We may admit that at a certain moment of time, which we can take as the present time to, 
the Hubble parameters along the orthogonal directions may be assumed to be approximately 
equal. Ha — Hb, even though we started with an anisotropic geometry. This hypothesis has 
been considered in for the case of a Kantowski-Sachs (KS) model. From this study was 
derived the conclusion that the classical tests of cosmology are not at present sufficiently 
accurate to distinguish between a FLRW model and the KS defined in that paper, with 
{HaQ — -fffto), except for small values of 6o- 

Taking Hag = H^^ for simplicity, one can then integrate equation ( P7| ) and find three 
different solutions, one for each k value. Figure 1 and Figure 2 shows the three kinds of 
behaviors as a result of integration. 

The behavior for the Kantowski-Sachs and Bianchi type-Ill cases depends on the ^Aq 
value. If VLao ^ ^Am' there will be a maximum value for y, (|/m), and since then ijiym) = 0, the 
slope of the curve x = x{y) will be infinite at that point. Specifically, we have i(?/m) = +oo 
for Kantowski-Sachs, and x{ym) = — cc for Bianchi type-Ill, even though x{ym) is finite. 
When ^Iam < ^Aq < ^Ac5 after x = y = 1 is reached we find an almost linear relation 
between the two scale factors x and y for the two models. While for Bianchi type-I model we 
have X = y ioT the present restrictions. So, we see that for the KS model, the scale factor a{t) 
starts from infinity if b{t) starts from zero. For the Bianchi type-I model, the scale factors 
are always proportional or even equal. In this situation we don't have an anisotropic model; 
in fact, we can easily prove that this model is isotropic by a properly reparametrization of 
the coordinates. For the Bianchi type-Ill model, the scale factor b{t) never starts from zero. 



"ao<"a- 




Figure 1: Scale factors relation, that is, the x y dependence for the three models Kantowski- 
Sachs (/c = 1), Bianchi type-I {k = 0) and Bianchi type-Ill [k = —1). We show the behaviour 
of x{y) when fl^^ < fiAM- Concretely we have for Kantowski-Sachs Qmo = 9 and ^Aq = 1.5; 
for Bainchi type-I Qmo = 2 and Q/y^ = —1; for Bainchi type-Ill Qmq = 1 and ^Aq = '~^■ 



^^<^A<^A. 




Figure 2: Scale factors relation, that is, the x y dependence for the three models Kantowski- 
Sachs (/c = 1), Bianchi type-I [k = 0) and Bianchi type-Ill [k = —1). We show the behaviour 
of x{y) when fiAj^j < ^a,, < ^Ac- The particular values for the plotting are, for Kantowski- 
Sachs Qmo = 2 and Qaq = 1-5; for Bainchi type-I Qmq = 0.5 and Q^,, = 0.5; for Bianch 
type-Ill Hmo = 0.2 and I^Aq = 0.4. 



'Ao 



but has an initial value different from zero when a is null. The following plot shows the zones 
in the 2-dimensional parameter space {fiMo,^Ao) where each model is allowed (Figure 3). 

Taking into account the analysis given in [Q, we may easily derive the qualitative be- 
havior oi y{t), since our equation (1l2|) is mathematically equivalent to his equation (3). Now, 
going back to Figure 1 and Figure 2, one can then determine the x{t). The plotting below 
summarizes the several possibilities for the three models: Kantowski-Sachs, Bianchi type-I 



a 




Figure 3: The Kantowski-Sachs model corresponds to the region above the straight hne 
{flko < 0); the Bianchi type-Ill model corresponds the region below the straight line; {Qko > 
0); the straight line represents the set of region for the Bianchi type-I model {^ko = 0). 

and Bianchi type-Ill models, respectively. 
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Figure 4: The scale factors x and y for the Kantowski-Sachs model [Q^o < 0) when ^Ao < 



^AJ^^. For the plotting we put fi 



Mo 



9 and ^2 



Ao 



1.5. 



The present technology allows us to "see" the epoch of last scattering of radiation at 

a redshift of about 1000, i.e., we can actually observe the most distant information that 

the Universe provides. The high level of isotropy observed from the Cosmic Microwave 

Background Radiation (CMBR) [|15| from this epoch to our days imposes that the two 

Hubble factors Ha and Hf, must remain approximately equal from this epoch to the present. 

In other words, we must impose a high isotropy level from the last scattering onwards, in 

our expressions, i.e., 

AH H„ - H, 



Hr, 



Ha 



^A<^A<^A. 




Figure 5: The scale factors x and y for the Kantowski-Sachs model [Qko < 0) when f^A^,^ < 
r^Ao < ^Ac- For the plotting we put Qmo = 2 and Q^o = 2. 



^a.<"a. 



° 




Figure 6: The scale factors x and y for the Bianchi type-I model {^ko = 0) when ^Aq < 0. 
For the plotting we put Qmo = 2 and Qa^ = — 1. 



such that 



AH 



^^ < 1. From COBE 4- Year data [|T6], 0, we have (fj ~ lO^^ and for last 

scattering epoch f-^j ~ 10^^. At the last scattering we may still consider Ha ^ H^ c::i H 
{H defined above). 

We computed several numerical integrations, with equation (|27|), by the following way: 
we gave values to ^2^0 and Qaq and integrated back in time, from now to the last scattering 
epoch. These Qmo and Qaq values were chosen such that, at last scattering epoch we had 



AH 



H 



1 - f^) = 1.7 X 10^^ or (^) = 1 ± 1.7 X 10"^ To make this we implemented 
a 8 order Runge-Kutta method ||TB|. 

We concluded that the sum Qmo + ^Aq must be close to the unity from above for 
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^A<^A<^A. 




Figure 7: The scale factors x and y for Bianchi type-I model {Qko = 0) when ^Aq > 0. For 
the plotting we put Qmo = 0.5 and fl\Q = 0.5. 
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Figure 8: The scale factors x and y for the Bianchi type-Ill model {Qko > 0) when Q^^ < 
f2AM- Fo^ the plotting we put Qmo = 1 and ^^Aq = — 1. 



Kantowski-Sachs and from below for Bianchi type-Ill models^. We summarized in the table 
below the result of imposing 



AH 



1.7 X 10 ^ for Kantowski-Sachs and Bianchi type-Ill 
models, supposing Hag = H^^^ (because (-^j ~ 10~^). 

From the table above we concluded that all combinations of Qq + Q^g near the unity 
are equally acceptable for reproducing a small anisotropy {{a/H)is ~ 10~^) at the last 
scattering. Nevertheless we paid special attention to the values of Qq ~ 0.3 and fl^g ~ 0.7, 
since they reproduce the better fit to recent observations |T9|. We have in this scenario 
\AH/H\ < 2 X 10^^ for Kantowski-Sachs and Bianchi type-Ill universes. All these models 



^It is obvious that for Bianchi type-I model {^Mo + ^Aq 
AH/ Ha = 0. 



1), with our restrictions, we have always 
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^A<^A<^A. 




Figure 9: The scale factors x and y for the Bianchi type-Ill model {Qko > 0) when ^Aa^ < 
r^Ao < ^Ac- For the plotting we put Qmq = 0.2 and ^Ao = 0.4. 



Table 1: Density parameters and relative difference between Ha and Hb for Kantowsky-Sachs 
and Bianchi type-Ill models. 



K-S 
K-S 
B III 
Bill 
B III 
B III 



n 



'Mo 



n 



'Ao 



J7o + ^Ao 



K-S 1 

K-S < 2 X 10' 



-15 



< 0.3 -h 7.0 X 10-9 

0.3 

1 - 10-10 

> 9.8 X 10-1^ 

> 0.3-10-11 
0.3 



< 1.7x 10"^ 
1 

0.7 

< 0.6 -F 7.0 X 10- 

>9.9 X 10-11 

1 - 10-13 

0.7 

> 0.7 -10-11 



1-F5.6 X 10-9 
1-F6.7X 10-1^ 
1 + 2.3 X 10-9 
1 + 2.3 X 10-9 
1 - 3.3 X 10-13 
1-6.7X 10-16 
1 - 3.3 X 10-12 
1-3.3 X 10-12 



^fco 






-1.7 X 10-*^ 
-2.0 X 10-1^ 
-7.0 X 10-9 
-7.0 X 10-9 
+1.0 X 10-12 
+2.0 X 10-15 
+1.0 X 10-11 
+ 1.0 X 10-11 



-1 
-1 
-1 
-1 

+1 
+1 

+1 

+ 1 



.6 X 10-*^ 
.4 X 10-6 
.7 X 10-6 
.7 X 10-6 
.3 X 10-6 
.3 X 10-6 
.8 X 10-6 
.8 X 10-6 



approach locally an exponential expanding state pO[ provided the cosmological constant if 



we consider i^A > ^Am- 



3 Conclusions 



For the Kantowski-Sachs model {VLkQ < 0) (see Figures 4 and 5), we conclude that if the scale 
factor h{t) starts from zero, then the scale factor a{t) will start from infinity and decreases 
afterwards. When f^Ao < ^\mi K't) reaches the maximum value recollapsing after that. 
So, a{t) will reach a relative maximum, when b{t) is maximum, (see Figure 4). After that, 
when b{t) = 0, ait) goes to infinity again. When ^Am < ^Ao < ^Ac? the scale factor b(t) 
grows indefinitely, giving place to an inflationary scenario. Then, a{t) decreases reaching a 
minimum value, and growing after that indefinitely, and becoming proportional to b{t) (see 
Figure 5). The initial singularity is of a "cigar" type. 

For the Bianchi type-I model {Qko = 0) (see Figures 6 and 7), the scale factors a(t) and 
b{t) are proportional or even equal. Thus, this model turns out to be an isotropic one (owing 
to our restrictions) and Qq + ^Ao = 1. However, when f^Ao < ^Aaj, ci{t) and b{t) reach 
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the maximum and recollapse after that. And when ^Am < ^Aq < ^Ac5 '^('^) ^^^ ^i^) grow 
indefinitely after an infiection. 

For the Bianchi type-Ill model {Qko > 0) (see Figures 8 and 9), when ^a^ < ^Am^ ^(t) 
starts from an initial non vanishing value {b{t = 0) = bo > 0), reaching a maximum and 
recoUapsing after that until reaches the same value for t = 0. Also, a{t) has a similar behavior, 
but starts from zero and recollapses to zero, nevertheless, a{t) exhibits a relative minimum 
when b{t) is maximum. When ^Am < ^Ao < ^Ac5 b{t) starts again from a non vanishing 
value (6o > 0), growing indefinitely with an infiection. In this case, a{t) starts from zero and 
grows indefinitely becoming approximately proportional to b{t). So, the initial singularity is 
of a "pancake" type. 

In conclusion, these models undergo isotropization becoming an asymptotically FLRW, 
except for the Kantowski-Sachs model (^2^0 < 0) with Qaq < ^Am ^^^ ^^ ^^e Bianchi type- 
Ill {flko >) with ^Aq < ^Am ■ Taking into account the accuracy of the measurements of 
anisotropy on one hand and the fact that we can always adjust the density parameters such 
that 1^0 + ^Ao ~ 1| = ^: with 6 ~ 10^^ on the other, we conclude that these models can still 
stand as good candidates to describe the observed Universe. 
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